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ABSTRACT 



Expendable Bathythermograph (XBT) temperature depth profile data 
from 17 locations in the north Pacific Ocean has been analyzed to extract 
information concerning the thermal microstructure. The microstructure 
information is then examined to determine if the distribution can be 
associated with the climatic regions of the north Pacific Ocean proposed 
by Tully [ 1964 ]. 

A modified binomial filter used by Roden [ 1971 ] to analyze STD data 
was applied to find a mean temperature series. A temperature per- 
turbation series results ^^hen the mean is subtracted from the original 
XBT sei'ies. Two resulting param.eters , average mean temperature 
perturbation vertical thickness scale and average maximum amplitude 
of the temperature perturbations , can be associated with the micro- 
structure. 

The average mean temperature perturbation vertical thickness scale 
was found to increase with latitude from 20°N to 55°N , and as the thick- 
ness scale increased, the maximum amplitude of the perturbations was 
found to decrease. 



2 



TABLE OF CONTENTS 



l. INTRODUCTION 10 

n. THEORY OF ACOUSTIC PROPAGATION IN A RANDOMLY 

INHOMOGENEOUS MEDIUM 13 

m. RECENT MEASUREMENTS OF MICROSTRUCTURE 24 

A. REVIEW 24 

B. MICROSTRUCTURE MECHANISMS 32 

1. Insolation 33 

2. Double Diffusion 33 

3. Internal Waves 34 

4. Current Shear 34 



IV. INFORMATION CONCERNING THE THERMAL MICROSTRUCTURE 



THAT CAN BE GATHERED FROM EXISTING EXPENDABLE 
BATHYTHERMOGRAPH (XBT) RECORDS 35 

A. AVAILABLE DATA 35 

B. EXPENDABLE BATHYTHERMOGRAPH SYSTEM 35 

C. DIGITIZING 37 

D. THERMAL MICROSTRUCTURE 41 

E. DATA SELECTED FOR ANALYSIS 47 

F. DESCRIPTION OF MICROSTRUCTURE ON XBT TRACE 48 



3 



G. OCEANOGRAPHIC REGIONS OF THE NORTH PACIFIC 

OCEAN 52 

1. Description 52 

2. Microstructure Distribution 55 

V. OPERATIONAL ANALYSIS 59 

A. OPERATIONAL REQUIREMENT 59 

B. OPERATIONAL METHOD 71 

C. POSSIBLE UTILIZATION OF MICROSTRUCTURE 

INFORMATION FROM AN XBT TRACE 72 

VI. CONCLUSION 73 

A. SUMMARY 73 

B. CONCLUSIONS 73 

C. RECOMMENDATIONS 74 

BIBLIOGRAPHY 76 

INITIAL DISTRIBUTION LIST 80 

FORM DD 1473 82 



4 







I 



I 



I 



I 



LIST OF TABLES 



I. Specifications of XBT System 36 

II. Subtropic Region 61 

ni. Subarctic Region 62 

IV. Transport Regions 63 

V. Other Regions 64 

VI. California Current Time Series 65 



5 



LIST OF FIGURES 



1. Coefficients of Variation of the amplitudes of direct signals 
plotted against the range to the receiver. Sound frequency, 

24 kHz [ after Sheehy, 1950 ]. 13 

2. Focusing Range (meters) as a Function of Frequency (kHz) 

and Patch Size (meters) [ Skudrzyk, 1957 ]. 18 

3. Increase in focusing due to spheriodal shape of a patch as a 

function of the ratio of the lengths of the vertical (b) to 
horizontal (a) axis [ Whitmarsh, Skudrzyk, and Urick, 1957 ].- 20 

4. Illustration of the various types of acoustic scattering due 

to a thermal patch in the sea [ Whitmarsh , Skudrzyk , and 
Urick, 1957 ]. 21 

5. Horizontal temperature distributions as measured in the sea 

[ after Whitmarsh, Skudrzyk, and Urick, 1957 ]. 22 

6. Temperature and temperature- gradient soundings for 0630 

24 August 1966 in the Mediterranean Sea [ Woods , 1968a ]. — 26 

7. Temperature and temperature-gradient soundings for 

Maltese waters, 14 September 1967 [ Woods, 1968b ]. 27 

8. Corrected data from the upper portion of the water column 

in the San Diego Trough on 18 September 1968. The curve 
on the right is the continuation of the curve on the left 
[Osborn, 1969 ]. 29 

9. Vertical profile of temperature under T-3 (84° 38'N, 

128° 21.6;'W, 19 March 1969). Shaded section of profile (a) 
shown as observed in profile (b) to left [ Neshyba , Neal , and 
Denner, 1969 ]. 30 

10. Time series taken from T-3 between 2145Z 17 March and 
2033Z 18 March 1969. Temperature increases from right 
to left [ Neshyba, Neal, and Denner, 1969 ]. 31 

11, . Unretouched Bathysonde profiles from the western Irish 
Sea of temperatui'e and conductivity showing pronounced 
microstructure [ Simpson, 1971 ]. 32 



6 



•■'■■vr 




12. Digitized record [Dale and Stevens , 1970]. 39 

13. Comparison of the original XBT trace with the trace 
reconstructed from the digitized data. The digitized 
trace is plotted to a scale of 90 meters/in due to 
computer requirements while the original trace is 

scaled to 92. 3 meters/in. 40 

14. Computer solution for modified binomial filter response 

with filtering intervals of 1,2, and 3 meters, 43 



15. Result of modified binomial filter applied to Theoretical 
XBT trace. From left to right. Theoretical XBT trace. 

Mean curve , and Temperature Perturbation curve for a 

7 meter wave length, 0. 5°C Temperature variation, and 
a 1 meter filtering interval. 45 

16. Result of modified binomial filter applied to theoretical 
XBT trace showing the amplification in the temperature 
perturbation amplitude when the wave length of the 
perturbation is on the slope part of the response curve 

of Figure 14. From left to right , Theoretical XBT trace , 

Mean cu 2 n/e , and Temperature Perturbation cur\'e for a 

34 meter wave length, 0. 5°C Temperature variation, 

and a 1 meter filtering interval. 46 

17. Composite diagram of the original XBT trace , the 
reconstructed Digitized trace , the Mean curve found by 
applying Roden's filter, and the Temperature Perturbation 



curve. The original trace is scaled to 92. 3 meters/in 

wdiile the remaining curves are scaled 90. 0 meters/in. 49 

18. Temperature perturbation curves for a filtering interval 

of 1, 2, and 3 meters. 51 

19. Regions of the North Pacific Ocean [Tully, 1964]. 53 

20. Drift in the North Pacific Ocean [Tully, 1964]. 53 

21. Plot of the Average Mean Temperature Perturbation 

Vertical Thickness Scale versus Latitude showing an 
increase in the vertical thickness scale for an increase 
in latitude between 20°N and 55°N. (A, points along 

Longitude 156-157°W; • , all other points). 66 



7 



22. Plot of the Average Maximum. Amplitude of the Temperature 
Perturbations versus Latitude showing the Amplitude de- 
creases with an increase in latitude between 20°N and 55°N. 

( A, points along Longitude 156-157°W; • , all other points). — 67 

23. Plot showing the Average Maximum Amplitude of the 

Temperature Perturbations decreases as the Average 
Mean Temperature Perturbation Vertical Thickness Scale 
increases between 20°N and 55°N. ( A , points along Longitude 

156-157°W; * , all other points). 68 



8 



ACKNOWLEDGEMENTS 



The author wishes to acknowledge with gratitude the guidance , 
patience, and understanding provided by Associate Professor W. W. Denner 
during the preparation of this thesis. 

The generous donation by Dr. N.E.J. Boston of his time and advice 
was essential for the orderly completion of this thesis. 

Special thanks are extended to Mr. Dean Dale and Mr. Paul Stevens 
of Fleet Numerical Weather Central, Monterey, California for their kind 
cooperation in providing data and assistance whenever requested. 

A final but most sincere word of appreciation is extended to Mr. 

J. F. T. Saur of the Fishory-Oceanogi'aphy Group of the United States 
Department of Commerce , National Oceanic and Atmospheric Administra- 
tion, National Marine Fisheries Service in La Jolla, California for providing 
copies of original XBT traces. 



9 



I. INTRODUCTION 



Temperature and salinity distributions in the ocean determine the 
propagation paths for acoustical energy. The gross temperature structure 
is used for ray path calculations to predict transmission characteristics 
[Cohen and Weinberg, 1969], However, sound measurements at one location 
reveal a fluctuating acoustic pressure level [Sheehy, 1950] attributed to 
variations in the fine temperature and salinity structure that cause small 
but important changes in the propagation paths. At higher frequencies 
the microstru.cture becomes increasingly significant. A theory’ for acoustic 
v.'ave transmission has been developed for a medium containing thermal 
microstructure assuming that the microstructure can be represented by 
inhomogeneities of a thermal nature that are spherical or spheroidal 
"patches" randomly distributed in space [For example, Liebemiann, 1951; 
Mintzer, 1953a, 1953b, 1954]. Tlie measured sound fluctuations have been 
related to the focusing or dcfocusing of acoustic waves by these "patches" 
at ranges less than ka^ (vdiere k is the acoustic wave number and a is the 
effective radius of the patch) and to interference at ranges greater than 
ka^ [Skudrzyk, 1957]. 

Recent findings concerning the microstructure do not support this 
"patch" distribution assumption. A layered microstructure with layers 
of small gradient a few to a few tens of meters thick separated by sheets 
of high gradient a few to a few tens of centimeters thick nd with 
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horizontal dimensions up to thousands of meters have been found in many 
regions of the oceans [Cooper, 1967; Cooper and Stommel, 1968; Neshyba, 
Neal, and Denner, 1969; Osborn, 1969; Roden, 1971; Simpson, 1971; 

Stommel and Fedorov, 1967; and Woods , 1968]. Several mechanisms have 
been proposed for generating the layered microstructure , and it may be 
that a different mechanism is dominant for different ocean regions. 

This layered microstructure concept requires a re-examination of the 
present theory for the transmission of acoustic energy in a randomly 
inhomogeneous medium. 

In 1957, VThitmarsh , Skudrzyk, and Urick concluded a paper on the 
scattering of sound in the sea with the following comment; ”. . . . the 
greatest g?.p seems to lie in the micrcthermal statistics of the ocean. 

If these were known better than they are now , apparently the amplitude 
statistics of the fluctuations of sound transmitted throu^ the sea could 
be predicted roughly, but with some assurance. Some 11 years later, 
during a meeting on Turbulence In The Ocean held at The University of 
British Columbia , R. W. Stewart "pointed out that a case existed for 
another degree of exploration of the ocean. While in the past, widespread 
surveys have sou^t to measure the temperature and salinity fields , we 
must now concentrate on a widespread search for fluctuations in the 

^ Whitmarsh, D.C. , Skudrzyk, E. , and Urick, R. J. , "Forward Scattering 
of Sound in the Sea and I ts Correl ati on with Temperature Microstruc- 
ture ," J_^_A£o^^_Soc^_Am. , V. 29, p. 1141, October 1957. 
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temperature and salinity fields. The only data that exists with such 
broad coverage is bathythermograph data. Magruder [1970] suggests that 
Expendable Bathythermograph (XBT) data might be useful. 

The purpose of this thesis is to see if information about the distri- 
bution of microstructure can be determined from XBT traces. If the 
microstructure distribution can be determined from an XBT trace , then 
using some of this available data , can the distribution be related to the 
climatic regions of the ocean for use as a rational basis for assessing 
the acoustic variability in these regions? Finally, what operational 
techniques can be developed using knowledge of the microstructure deter- 
mined from the XBT trace that would be of benefit to the Navy in the 
operating environment? 

Tens of thousands of XBT traces are on file at the National Oceano- 
graphic Data Center (NODC) , with the number increasing each day. Under 
the Fleet Numerical Weather Central (FNWC) , Cooperative Oceanographic 
Observation Program , XBT traces are also available from commercial 
carrier trans-atlantic and trans-pacific routes. While the distribution 
is by no means uniform in space or time , it is the only source of micro- 
structure data covering the world oceans. 



^ Osborn, T.R. , Report of the Symposium on Turbulence in the Ocean , 
Institute of Oceanography, University of British Columbia, Vancouver, 
British Columbia, p. 20, July 1970. 
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II. THEORY OF ACOUSTIC PROPAGATION IN A 
RANDOMLY INHOMOGENEOUS MEDIUM 



Sheehy [1950] pubUshed data showing that the magnitude of the 
acoustic pressure maximum amplitude fluctuations (Coefficient of 
Variation) for 24 kHz direct path transmissions increased proportional 
to the square root of the range. However, it should be noted from 
Figure 1 that the scatter of data is rather large. 




RANGE, YARDS 



Figure 1. Coefficients of Variation of the amplitudes of direct 
signals plotted against the range to the receiver. Sound frequency, 
24 kHz [after Sheehy, 1950]. 
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Bergmann [1946] using the application o£ ray optics to investigate 
the propagation o£ radiation in a medium with random inhomogeneities , 
derived an approximate formula which showed that the signal level fluctua- 
tions increased in proportion to the three^halves power of the range. He 
pointed out the variance with Sheehy's results , but offered no explanation 
for the differences. 

Liebermann [1951] , using an autocorrelation function calculated from 
measured horizontal small-scale temperature inhomogeneities , showed 
that the refractive effects of the horizontal inhomogeneities produce a 
local warping (focusing or defocusing of sound) that causes the acoustic 
intensity to fluctuate. He found that refraction should only occur if the 
frequency is high enough to make the dimensions of the inhom.ogeneities 
greatly exceed the wave length of the radiation. 

Mintzer [1953a] derived a first order approxim.ation for the acoustic 
pressure amplitude variation of an acoustic pulse traveling in a medium 
with small random variations in the refractive index that are a result of 
random temperature variations. Assuming that the pressure amplitude 
fluctuations are small enough to allow the use of a single scattering 
approximation and that the refractive index at a point only changes during 
the time between pulses and not while the pulse is passing the point, he 
solved the Born approximation to the wave equation for the pressure 
amplitude of the acoustic pulse. Then using the Coefficient of Variation, 
V , defined as : 
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v^= 

l<p>l ^ 



where 



<> = ensemble average 
P = mean presure amplitude 

and by averaging the pressure amplitude over the pulse length instead o£ 
using the maximum value o£ the pressure amplitude as Sheehy did with his 
experimental data, Mintzer’s theoretical analysis yields the £ollowing. 
£or ranges large compared with the correlation distance: 

V^= 2ko a^r p)dp , where 

*o 

N( p)= autocorrelation £unction o£ the retraction 
index variation taken along the line between 
source and receiver 
k(j= wave number o£ the sound 
r = range £rom source to receiver 
p= rms retractive index variation 
Using a Gaussian correlation tunction the coetticient ot variation 
becomes: 

J- 2 2 % % 

V = ( TT^ k^aa) r , where a = correlation 
length 
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This theoretical result is in agreement with Sheehy's experimental 
results [Figure 1], 

Later Mintzer [1953b] determined over what range the first assump- 
tion (use of a single scatterer approximation) was valid. He found that 
for large ranges compared with the acoustic wave length and the inhomo- 
geneity size , and a short pulse length compared to the time required for 
a refractive index change: 

(1) The Coefficient of Variation of the pressure 
amplitude is: 

^ 1<Q O' ^ ar , ka » 1 
= 7 t'^ ko a ^ a" r , ka « 1 

(2) The single-scattering assumption used is valid 
for: 

ko r « (k^a a ^ 

Mintzer also made a comparison with ray theory and found that for a 
frequency high enough to satisfy: 

rA<-Q3.^« 1 

the intensity coefficient of variation computed by wave theory and ray 
theory are the same. 
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In a third paper, Mintzer [1954] considered the second basic assump- 
tion, that the refraction index varies only during the time between pulses. 
He found that "if the pulse length is much less than the time constant of 
the refractive index correlation function , the coefficient of variation is 
independent of the pulse length , and the pulse -correlation function is 
found to be identical with the time-dependent part of the refractive index 

O 

correlation function. 

Skudrzyk [ 1957] discussed the pressure amplitude fluctuations caused 
by water inhomogeneities considered to be spherical "patches" that are 
always large in comparison with the acoustic wavelength. He showed that 
in the vicinity of the receiver pressure amplitude fluctuations are caused 
by focusing or defocusing from the lens action of the patches , and beyond 
a certain range the fluctuations are caused only by an interference field 
which results from deflected rays arriving at a point with phase shifts 
due to unequal paths of travel. He also pointed out that the interference 
component of scattering is independent of the shape of the patches , while 
the focusing component can be significantly affected by the shape. For 
example , a spheroidal patch having a vertical axis less than the horizontal 
axis will have a larger focusing component than a spherical patch. At low 
frequencies (below 20 kHz) interference scattering dominates , while at 

Mintzer, D. , "Wave Propagation in a Randomly Inhomogeneous 
Medium. Ill," T. Acous. Soc. Am. , v. 26, p. 186 , March 1954. 
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higher frequencies the interference component may be so small in 
comparison with the focusing component that at ranges greater than the 
focusing range the Coefficient of Variation is almost constant. The 
focusing range ( r^ ) is the range beyond wliich the focusing does not 
change and additional pressure amplitude fluctuations are caused by 
interference only. Skudrzyk found the focusing range to be a function 
of frequency and path size. 

k = acoustic wave number 

•y 

r^ = k a"^ , where 

a = effective radius of patch 

This result is pictured graphically in Figure 2. If the sonar frequency 

and the patch size are known , then the focusing range can be determined. 

(For example, a 24 kHz sonar has a focusing range of 200 meters, if the 

patch size is assumed to be 140 centimeters). 

a ► 




Figure 2. Focusing Range (meters) as a Functii ' of Frequency 
(kHz) and Patch Size (meters) [Skudrzyk, 1957]. 
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Skudrzyk shows that the Coefficient of Variation of the pressure 
amplitude due to the focusing effect, V£, is; 



V 






where « is the relative root-mean- 
square deviation of sound velocity 



while the Coefficient of Variation due to the interference effect is; 




, 2 2 

2 ^ a a (r-r^ ) 



Beyond the focusing range , scattering is due to both focusing and inter- 
ference so the Coefficient of Variation becomes; 

v^ =vl+ 

p t 1 

and for spherical patches ; 

V^= k^a (r-0.44ro) 

which is Mintzer's [1953b] result except for a factor of 1/2. If the 
patches are spheroidal rather than spherical , then the focusing may be 
increased significantly. Whitmarsh, Skudrzyk, and Urick [1957] found 
the focusing to increase by a factor 

4 6 

Vp (Spheroidal) ^ + 6 

Vp (Spherical) 16 

X is the ratio of vertical to horizontal axis of the 
patch 
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Figure 3 depicts this ratio graphically. The interference component of 
scattering is virtually independent of the vertical axis dimension. 




Figure 3. Increase in focusing due to spheroidal shape of a 
patch as a function of the ratio of the lengths of the vertical 
(b) to horizontal (a) axis [Whitmarsh, Skudrzyk, and Urick, 
1957]. 



Since published experimental results pertained only to one frequency 
and did not include temperature and acoustical fluctuation measurements 
made as one observation, Whitmarsh, Skudrzyk, and Urick [1957] conducted 
a series of experiments using 25 kHz and 60 kHz transmissions to provide 
additional data to test the existing theories for scattering of acoustical 
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energy in the sea. They used a gradiometer bar suspended fronri a buoy 
(a 12 foot long bar with 9 thermistor beads) to record temperature 
gradients as small as 0. 001°C and to record the temperature fluctuations 
as a function of time. They point out that all the existing theories for 
scattering of sound in an inhomogeneous medium are based on the assump- 
tion of a spherical patch thermal inhomogoieity (cane such patch is illustra- 
ted in Figure 4). The theories further assume that the patch distribution 
in the medium can be represented by an exponential or Gaussian correlation 
function. 




Figure 4. Illustration of the various types of acoustic scattering 
due to a thermal patch in the sea [Whitmarsh , Skudrzyk , and Urick , 
1957]. 



The teinperatiBre microstructure measurements revealed that the assumed 
correlation function is invalid and in fact the results she the temperature 
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fluctuations increase vvath the 2/3 power as the thei-mistor spacing 
increases. This is the prediction of the Kolmogorov theory for homo- 
geneous turbulence, see Figure 5. However, even replacing the correlation 




Figure 5. Horizontal temperature distributions as measured in the 
sea [after Whitmarsh, Skudrzyk, and Urick, 1957]. 



function by a new function representative of the Kolmogorov distribution, 
the theories for sound transmission can still be considered vaUd. 
Whitmarsh , Skudrzyk , and Urick concluded from their experiments that 
the present theories are capable of predicting the approximate variation 
in the pressure ampHtude to be expected for acoustical transmissions in 
an inhomogeneous medium. 
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Sagar [1960] analyzing the experimental data he collected o££ New 

Zealand suggests that sound intensity fluctuations which obey the Mintzer 

2 2 2 

dependence on range oc r£ a a) are related to insolation a££ects. 

His data is interpreted to show a diurnal cycle o£ grcwth and decay o£ 
tile temperature microstructure in the surface layers o£ the sea with 
the structure most developed by late afternoon and then disappearing over- 
night. He also found and described the water to have a "sandwich” or 
layered structure rather than the classical "patch" structure. 
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m. RECENT MEASUREMENTS OF MICROSTRUCTURE 



A. REVIEW 

Stommel and Fedorov [1967] used a salinity, temperature, and depth 
recorder (STD) to take measurements o£ the fine structure (nucro- 
structure) in the main thermocline off Timor in the Arafura Sea and 
Mindanao in the Philippine Sea. They were investigating reports of fine 
structure which began to acciunulate with the advent of an instrument 
like the STD, which continuously records temperature and salinity as a 
function of depth. They found the vertical microstructure in the main 
thermocline to consist of hundreds of layers (].aminae) of rather homo- 
geneous temperature from 2 to 40 meters thick separated by a thin layer 
(2 to 5 meters thick) of high temperature gradient. As an example , at 
one station a mean thermal gradient of 0. 040°C/meter was found, while 
in a homogeneous layer the gradient was greater than 0. 35°C/meter. By 
tracing some of the well defined laminae in a horizontal direction , they 
found the laminae extended from 2 to 20 kilometers. Though the STD 
resolution was only 2 meters, Stommel and Fedorov felt the micro- 
structure existed on a smaller scale and that further investigation would 
be of interest. 

Cooper and Stommel [1968] made microstructure measurements in the 
main thermocline off Bermuda using an STD. They too found a layered 
structure with the homogeneous layers from 3 to 5 meters thick and the 
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gradient layers between 10 and 15 meters thick. In the layers of 
hi^ gradient, temperature changes as high as 0. 50°C were found. The 
horizontal extent of the layers was between 400 and 1000 meters. This 
type structure was present out to a distance of 120 kilometers from 
Bermuda. Based on their measurements , Cooper and Stommel suggested 
that most of the Sargasso Sea exhibited the same characteristics in the 
main thermocline. 

Woods [1968a] conducted an extensive study of the microstructure 
in the summer thermocline of the Mediterranean Sea. A temperature 
gradient meter was used to obtain measurements of temperature gradient 
while being lowered and temperature while being recovered. The instru- 
ment consisted of two thermistors mounted 50 centimeters apart on a 
vertical fiber glass rod. Only one of the thermistors was used for the 
temperature measurements. The instrument had an 0.01°C accuracy and 
a depth limitation of 50 meters. Woods found layers with a moderate 
temperature gradient that were thick in comparison with alternate layers 
(sheets) of high temperature gradient. Figure 6 is a plot of his results 
for one sounding. High gradient sheets can be seen at depths of 10 , 19 , 
24, 32, 41, and 50 meters. By injecting dye from suspended packets into 
the different layers , Woods was able to associate laminar flow to thin 
high gradient layers and a weak turbulent flow to the thick homogeneous 
layers. The results of these two different type flows is a strong sheer 
across the hi^ gradient sheets. 
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Figure 6. Temperature and temperature- gradient sounding 
for 0630 24 August 1966 in the Mediterranean Sea [Woods, 1968a], 

The next year Woods [1968b] conducted another investigation of 
vertical heat flow through the upper ocean in the Mediterranean Sea off 
Malta. Again the sheets in the thermocline were found as can be seen 
in Figure 7. 
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Figure 7. Temperature and temperature-gradient soundings 
for Maltese waters, 14 September 1967 [Woods, 1968b], 

Lovett [1968] took temperature profiles in the California current 
off Baja California using a sound, velocity, temperature, and pressure 
recorder (SVTP). He too found a definite layered structure in all the 
temperature profiles. He suggests that the sharp gradients found in 
the step- like structure may have an effect on acoustic scattering, 
especially so at high frequencies. 
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Osbom [1969] using a free fall instrument called a thermal 
gradiometer made microstructure measurements in the San Diego Trou^ , 
Florida current and in the central Pacific between San Diego and Hawaii. 
He found regions of high gradient to be about 1 meter thick and have a 
gradient around 10 times the mean gradient. Figure 8 is a record of the 
results he obtained with the gradiometer. By tracing a region of high 
gradient horizontally he found they extended up to 750 meters if the 
region is in the seasonal thermocline , but only a few hundred meters if 
the region is at a depth greater than 400 meters. Osbom suggests that 
the micros tmcture is on a scale of centimeter order. 

Neshyba, Neal, and Denner [1969] made measurements of the temper- 
ature mirr-ostruclure in the Arctic Ocean using an XBT lowered mechan- 
ically from Fletcher's Ice Island (T-3). The microstructure found had 
a step-like structure with layers that were extremely homogeneous as 
can be seen in Figure 9. The transistion layer was found to be less than 
20 centimeters thick. Figure 9b shows a section of the trace recorded 
at high gain (the lip at the end of a transition layer is due to recorder 
overshoot). The steps from 220 to 300 meters were found to be very 
stable. Figure 10 is a time series taken with the XBT modified to fall 
at 7 ft/sec instead of the normal 20 ft/sec, and demonstrates the stabil- 
ity. Below 350 meters , the homogeneous layers are larger, but still are 
separated by steps of temperature change. This is also apparent from 
Figure 10. It should be noted that these XBT recordings are independent 
of ship motion since the ice island is a very stable platform. 
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Figure 8. Corrected data from the upper portion of the 
water column in the San Diego Trou^i on 18 September 1968. 
The curve on the right is the continuation of the curve on 
the left [Osborn, 1969]. 
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TEMPERATURE (degrees Celsius) 




Figure 9. Vertical profile of temperature under T-3 
84 38'N, 128 21. 6' W, 19 March 1969). Shaded section of 
profile (a) shown as observed in profile (b) to left [Neshyba, 
Neal, and Denner, 1969]. 
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DEPTH (meters) 




\ 2145 ( 3 / 17 / 69 ) 0300 ( 3 / 18 ) 



2033 ( 3 / 18 ) 



Figure 10. Time series taken from T-3 between 2145Z 
17 March and 2033Z 18 March 1969. Temperature increases 
from right to left [Neshyba, Neal, and Denner, 1969]. 
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Simpson [1971] , using a Bathysonde CTD system , examined the density 
and thermal microstructure in the western Irish Sea. He found vertical 
layers up to 10 meters thick and in cases , less than 1 meter thick. The 
average gradient was as low as 0. 002°C/meter in layers up to 10 meters 
thick and as high as 2*^C/meter in the layers of sharp gradient. Figure 11 
gives three temperature profiles of the region. Horizontally the layers 
extended up to a distance of 3 kilometers. 




Figure 11. Unretouched Bathsonde profiles from the western 
Irish Sea of temperature and conductivity showing pronounced 
micro structure [Simpson, 1971]. 



B. MICROSTRUCTURE MECHANISMS 

The exact nature of the processes producing, maintaining, or destroy- 
ing the thermal microstructure still are unclear because of the complexity 
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of the problem. Some mechanisms have been proposed to explain the re- 
sults of observations and experiments. 

1. Insolation 

Insolation (radiation from the sun and atmosphere) results in the 
heating of the water's surface. The heat is transported downward mainly 
by mixing caused by wave action. The wnd-wave action tends to eliminate 
the surface layer microstructure. However, Sagar [1960] suggests that 
after sunrise the insolation may be large enough to dominate over wave 
action, and so a layered structure could be transported down the thermal 
gradient. Tlie depth that the layers reach is dependent on wave action. 
Waves are required to produce the mixing, but if they are too strong, 
then vigorous mixing will eliminate the layered structure. Sagar [1960] 
found that if cloud cover existed or if it was during a period of darkness , 
the microstiucture was absent. 

2. Double Diffusiot) 

This mechanism results from the differences in the molecular 
diffusivity of temperature and salinity which are related by an approxi- 
mate 100:1 ratio. Due to the differing diffusivity in regions where a 
warm saline mass of water overlies a cooler mass , the masses become 
unstable, forming a layered structure. Much evidence from ocean obser- 
vations [Cooper and Stommel , 1968 , and Tait and Howe , 1968] , a theor- 
etical study [Turner, 1968], and laboratory experiments [Turner and 
Stommel, 1964, Turner, 1967, and Stem, 1969] support this mechanism. 
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Stem [1969] feels that a turbulent mixiRg process is taking place in thick 
layers that are separated by thin discontinuity layers. 

3. Internal Waves 

Woods [1968a] observed internal waves along an individual micro- 
structure sheet to be in long coherent trains. If these waves become 
unstable , they would break , giving up their energy to turbulent processes 
causing the sheets to be eliminated by mixing. Orlanski and Bryan [1969] 
have proposed that internal waves breaking could also establish a step struc- 
ture by leaving a ''blob" of water behind that would spread out into a thin 
layer. 

4. Current Shear 

Current shear can cause vertical turbulent mixing when the 
Richardson Number is less than 0. 25 [Miles and Howard, 1963]. This 
turbulence occurs as a patch which is smeared out , leaving a scar that 
can lead to microstructure. 

Lovett [1969] and Pingre [1969] have proposed that these isolated 
turbulent patches can also be spread horizontally by turbulence and thus 
form a sheet. Dye studies illustrate this layering in the deep ocean 
[Schuert, 1969]. 

The complexity of the problem is apparent when it is realized 
that not one , but probably a combination of these mechanisms is occurring 
simultaneously. 
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IV. INFORMATION CONCERNING THE THERMAL MICROSTRUCTURE 
THAT CAN BE GATHERED FROM EXISTING EXPENDABLE 
BATHYTHERMOGRAPH (XBT) RECORDS 

A. AVAILABLE DATA 

Introduction o£ the XBT system provided a method relatively 
inexpensive , yet at the same time accurate , o£ obtaining temperature- 
depth pro£iles originally developed £or ASW operations , the system's 
application has been expanded to include Navy research ships and merchant 
ships equipped under the Fleet Numerical Weather Central (FNWC) , 
Cooperative Oceanographic Observation Program. The commercial ship 
routes o£ the Paci£ic and Atlantic Oceans are a principal source o£ ocean 
temperature structure data. XBT's taken during a cruise are £orwarded 
at the end o£ the cruise to FNWC , Monterey £or digitizing and then 
archived £or JBurther use. One copy o£ the digitized record along with the 
original XBT traces is sent to the National Oceanographic Data Center 
£or permanent storage in their archives. As o£ May 1971, approximately 
60,000 XBT's had been digitized by FNWC , Monterey and some 28 ,000 o£ 
these were £rom the northern Pacific area [Dale, 1971]. 

B. EXPENDABLE BATHYTHERMOGRAPH SYSTEM 

The expendable bathythermograph system currently in use is manu- 
factured by the Sippican Corporation of Marion, Massachusetts. The 
system is composed of the probe (the expendable part of the system) , 
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a launcher, and a recorder. It is an automatic system for measuring 
and recording sea temperature as a function of depth. The specifications 
for the system are given in Table I [Sippican Corporation, 1968]: 

TABLE I 



Temperature range 

Temperature accuracy 

Depth range 

Depth accuracy 

Deployment cycle time 
Gradient resolution 

Full scale stylus travel 



28°F to 96°F 
-2°C to 350 c 

+ 0. 4°F 
+ 0. 2°C 

0 to 1500 feet 

0 to 2500 feet 

+ 2% or 15 feet , whichever 
is greater 

90 seconds for 1500 ft drop 

63% of a step change in 
temperature in 3 feet; 

95% of a step change in 
9 feet 

1 second 



The probe is released from a cannister which pays out wire to 
compensate for the movement of the launcher platform while the probe 
itself pays out wire to compensate for its fall. The probe falls with a 
velocity of 20 ft/sec , and has been found to have an average terminal 
velocity of 19.9 +0.7 ft/sec [Gouzie, Sanders, and Littlehale, 1966]. 

A thermistor is mounted in the nose of the probe , and as the probe 
descends through the water column, the water temperature is sensed 
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as it flows past the thermistor and out the tail of the probe. The 
response time for the thermistor is 110 milliseconds nominal and 130 
milliseconds maximum [Choate, 1970]. The temperature resolution 
of the probe is dependent on the magnitude of the temperature change 
over a certain depth interval. 

The recorder uses a pressure sensitive paper for recording the 
temperature profile. The paper moves vertrically at a constant rate 
proportional to the fall rate of the probe. The chart paper displays 
311. 656 ft/in or 92. 307m/in and 9. 714°F/in or 5. 286°C/in, giving a 
resolution of 2 ft and 0. 1°F for visual redding. 

C. DIGITIZING 

The digitizing process used by FNWC, Monterey is completely 
described by Bauer [1969]. The CALMA Model 480 digitizer is a X - Y 
recorder that records at least every 0.15 inch of movement in either 
the X or Y direction. The speed the stylus is moved down the trace 
determines the number of points recorded. The maximum resolution 
is 0. 01 inches of movement in the X or Y direction which is equivalent 
to . 053°C or . 097°F in temperature and 0. 92 meters or 3. 12 feet in 
depth. Each trace is digitized twice. 

The digitized data is further processed using a DIGITBT computer 
program. The program first puts the data in a common format.- Then 
the two digitized records are compared point by point in the X and Y 
dimensions. If the maximum differences in both directions are less 
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than 0. 035 inches (0. 186°C or 0. 340°F in temperature and 3. 22 meters 
or 10. 92 feet in depth) , the two records are averaged on a point by point 
basis to produce an average trace which is used for the remainder of the 
DIGITBT program. If the deviation in either direction is greater than 
0. 035 inches , then the records are rejected and must be redigitized. 

The program next selects points that will be printed out and stored 
on the XBT record tapes. These points are called INFLECTION POINTS. 
Both the first and last points of the record are automatically selected 
as inflection points. The other points of the record are screened in 
’Tjatches". For each batch the observation is found that is most differ- 
ent from that expected by linear interpolation between the end points 
of the batch. If that difference exceeds 0.015 inches (0. 079°C or 
0. 146°F and 1. 38 meters or 4. 68 feet) then the observation is selected 
for output. Otherwise the batch expands to include the next value and 
the inflection test is repeated until an inflection is found and until the 
entire record has been screened. Any point that can be reconstructed 
by linear interpolation is not required and therefore is removed from 
the output. The DIGITBT program converts the inflection points from 
inches to meters and Centigrade degrees for output. The depth is 
rounded to the nearest whole unit and temperature to the nearest tenth 
of a degree. Figure 12 is an example output from the DIGITBT program 
and is the information stored in FNWC , Monterey's records. The title 
line contains a ship identification code , day , month , year , time , 
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Figure 12, Digitized record [Dale and Stevens, 1970 ], 



latitude , longitude , base line temperature , grid type , and maximimt 
depth o£ probe (15 for 1500 ft XBT or 25 for 2500 ft XBT). Figure 13 
is a comparison of the original XBT trace and the reconstructed trace 
from, a DIGITBT record sim,ilar to the exam.ple in Figure 12. As co'old 
be expected, the digitizing process has eliminated some of the small 
scale microstructure. 

For this thesis the data used is from the FNWC, Monterey's data 
bank that has been processed using the DIGITBT program. 
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Figure 13. Comparison o£ an original XBT trace with the trace 
reconstructed from the digitized data. The digitized trace is plotted 
to a scale of 90 meters/in due to computer requirements \^^ile the 
original trace is scaled to 92. 3 meters/in. 
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D. THERMAi MICROSTRUCTURE 

Thermal microstructure for the purpose of this thesis is defined 
as the temperature fluctuations as a function of depth that exist on a 
scale that is smaller than the scale of the gross temperature structure 
which is derived from the large scale processes such as weather , latitude , 
season, and water mass transport. 

To analyze the XBT trace for microstnicture , a method was sougjit 
to separate the temperature perturbations (microstructure) from the 
mean of the XBT trace. A method to filter the trace and subtract the 
mean portion of the curve was desired. Several type filters were consid- 
ered. A binomial filter was selected primarily because Roden [1971] had 
applied a modified binomial filter to similar data with favorable results. 

Roden utilized a low pass filter to develop a mean series. This mean 
series was then subtracted from the original series to give a temperature 
perturbation series. The infrmation was therefore effectively passed 
through a high pass filter. Using this scheme the temperature pertur- 
bations are considered to be associated with the thermal micro structure. 
If the original series is denoted by s(z), then the mean series, s*(z), is 
fovind by: 

m 

2 

s*(z)=^ w(i)s(z+l) 

1 — -Y 
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where w(i) is a smoothing array that has the following properties: 

w(i ) = w(-i) 

^w(i) = 1 

i 

The smoothing array is found by evaluating the below modified binomial 
sequence. 




The parameter, m, is even and controls the width (number of weighted 
points) of the array. The parameter, p, is an arbitrary exponent which 
controls the peakedness of the array. Roden found that a statistically 
stable mean series could be obtained if m = 40, with p = 0. 5. Increasing 
m produces more smoothing, but the series is physically shortened as will 
be seen later. The response of the modified filter, r(k), is given by 

r(k) = l-^w(i)cos(ik d S ) 
i 

where k = wave number ( "^Vwave length) and = filtering interval 
(distance between the 41 weighted points). Figure 14 is the solution for 
the response function for filtering intervals ( ) of 1, 2, and 3 meters. 
The smaller A the sharper is the response curve at a cut-off wave 
length, however, the cut-off wave length is lowered. It should be 
remembered that with a binomial filter the amplitude response is not very 
predictable in the region of cut-off. 
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Figure 14. Computer solution for modified binomial filter 
with filtering intervals of 1, 2, and 3 meters. 
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A theoretical XBT trace was constructed for test purposes. An 
isothermal layer of 50 meters was followed by a thermocline simulated 
with an exponential decay function (e that decayed in 650 meters. 
The microstructure was constructed by superimposing a sine wave 
function (sin2 t kz) on the exponential function. The amplitude was 
varied to represent the temperature fluctuations desired, and the wave 
length varied to represent the microstructure layer thickness scale. 
Figure 15 is a composite picture of a theoretical trace , the mean 
filtered trace , and the resulting temperature perturbation curve for a 
7 meter wave length, 0. 5°C temperature variation, and a filtering 
interval of 1 meter. As expected , the perturbation curve is the input 
sine wave curve with both the wave length and amplitude of the temper- 
ature variation being preserved. Figure 16 is another set of theoretical 
curves showing that for the one meter filtering interval if the wave 
length is 34 meters and the temperature amplitude variation is again 
0. 5°C the temperature perturbation amplitude has been ampEfied, while 
the wave length is preserved. It should be noted that the 34 meter wave 
length is on the slope part of the response curve in Figure 14, and this 
shows the instability of the filter in this region as predicted. If the 
filtering interval is increased to 3 meters then the original sine wave 
amplitude will again be passed. 
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Figure 15. Result o£ modified binomial filter applied to Tbeorectical 
XBT trace. From left to right. Theoretical XBT trace, Mean curve, 
and Temperature Perturbation curve for a 7 meter wave length , 0. 5°C 
Temperature variation, and a 1 meter filtering interval. 
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Figure 16. Result o£ modified binomial filter applied to theoretical 
XBT trace showing the amplification in the temperature perturbation 
amplitude when the wave length of the perturbation is on the slope 
part of the response curve of Figure 14. From left to right, Theo- 
retical XBT trace , Mean curve , and Temperature Perturbation curve 
for a 34 meter wave length, 0. 5°C Temperature variation, and a 
1 meter filtering interval. 
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E. DATA SELECTED FOR ANALYSIS 



A selection o£ 129 XBT traces has been analyzed for this thesis using 
Roden's modified binomial filter with 1,2, and 3 meter digitizing inter- 
vals. All of the traces were from the northern Pacific Ocean. Digitized 
information was furnished by the Fleet Numerical Weather Central (FNWC) , 
Monterey , and was processed using the IBM 360/67 at the Naval Post- 
graduate School's W. R. Church Computer Center. The traces were 
extracted from the FNWC data bank by selecting geographic grid squares. 
The selection of these grid squares was based on three criteria. The first 
and most important was to get a sampling from the Subtropic , Subarctic , 
Kuroshio Transport and Region of Confluence climatic regions proposed by 
Tully [1964]. Secondly it was desired to secure data for comparison with 
copies of original XBT traces obtained from J. F. T. Saur of NCAA's 
National Marine Fisheries Service at the Fishery-Oceanography Center in 
La Jolla, California. The third criteria was to request a reasonable quan- 
tity of data from FNWC. Data was desired for all months of the year in 
each grid square , but unfortunately the coverage fell short of expectations. 
Another shortcoming in the data is that the Subarctic region has not been 
sampled as often with XBT as the Subtropic region and therefore , data is 
not very plentiful from this area. 
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F. DESCRIPTION OF MICROSTRUCTURE ON XBT TRACE 



Figure 17 is a composite diagram o£ an actual XBT trace, the 
digitized trace, a mean curve found with Roden’s filter, and the temper- 
ature perturbation curve using a 1 meter filtering interval. The scale of 
the actual XBT trace is 92. 3 meters/in while the other three traces are 
scaled 90 meters/in for computer requirements. By comparing the actual 
XBT trace and the digitized reproduction it can be seen that some of the 
very fine scale temperature microstructure has been eliminated by the 
digitizing process and some larger scales appear to have been amplified 
by the process. The exactness of the digitized trace is considered to be 
sufficiently reliable to achieve results that are realistic and comparable. 

The first and last 20 meters of the mean trace are lost because the 
filter used 20 points zt 1 meter intervals on both sides of a depth to 
compute the mean value for that depth. If a 2 meter filtering interval 
is used, then a total of 80 meters is lost and the 3 meter interval reduces 
the length by 120 meters. Since the XBT record is comparatively short 
to begin with, a 750 meter trace would be reduced by 16 percent when using 
the 3 meter interval, and if the trace is a 1500 ft XBT (460 meters) 26 
percent of the record will be lost. 

By closely comparing the digitized XBT trace and the temperature 
perturbation curve in Figure 17 , it can be seen that any time an inflection 
point (as defined in Section IV C) is encountered on the XBT trace , a change 
in the slope results on the temperature perturbation curve. The number of 
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Figure 17. Composite diagram o£ the original XBT trace , the 
reconstructed Digitized trace , the Mean curve found by applying 
Roden’s filter, and the Temperature Perturbation curve. The 
original trace is scaled 92. 3 meters/in while the remaining cu^'es 
are scaled 90.0 meters/in. 
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slope changes can be counted and di\dded into the length o£ the trace 
producing a mean temperature perturbation vertical thickness scale which 
is representative o£ the microstructure mean vertical scale. The same 
value can be obtained by dividing by the number o£ in£lection points that 
was selected £or the digitized record. To be sure that only the micro- 
structure is being considered, the inflection points or slope changes 
above the seasonal thermocline must not be used nor the depth interval 
from the surface to the seasonal thermocline. This must be determined 
for each trace independently. The amplitude of the temperature pertur- 
bations is defined as the magnitude of the temperature difference between 
the mean curve and the XBT trace. 

Figure 18 demonstrates how the t emperature perturbations are 
affected by changing the filtering interval. Examination of the three 
traces below 100 meters shows that if the filtering interval is increased, 
then the shape of the perturbation curve is changed. In general , it was 
found that the mean temperature perturbation vertical thickness scale 
was increased when the filtering interval ( A$) was increased, and the 
temperature amplitude was increased, but not at all depths. It should 
also be noted that more of the trace is influenced by the mixed layer depth 
as the interval is increased, so that for the 3 meter curve 150 meters 
are not useable in computing the mean temperature perturbation vertical 
thickness scale. 
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Figure 18. Temperature perturbation curves for a filtering 
interval of 1 , 2 , and 3 meters. 
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Previous observations of the microstructure indicate scales below a 
few meters are important. Since increasing the filtering interval pro- 
duced an increase in the mean temperature perturbation vertical thickness 
scale , the perturbation amplitude and the variance , it was decided that 
the 1 meter filtering interval was more effective for processing the 
digitized XBT data. 

G. OCEANOGRAPHIC REGIONS OF THE NORTH PACIFIC OCEAN 
1. Description 

As previously mentioned , Tully [1964] divided the north Pacific 
Ocean into seven regions to be used for examining and forecasting oceano- 
graphic characteristics and behavior on a similar scale as the Weather 
Services have been doing for many years. In defining the regions , he 
considered volume transport, heating and cooling, precipitation and evap- 
oration, mixing (wind and convection), internal waves, and the effects 
these processes have on changing the properties of sea water. Figure 19 
is a diagram of the Pacific Ocean and the seven regions. The Tropic, 
Subtropic , Arctic , and Subarctic regions are the water mass formation 
areas where due to slow and/or zonal transport, the water mass has time 
to establish the regional properties. The normal time required to achieve 
these properties is one heating and cooling season. 

The Tropic region is bounded by Latitudes 10°S and 10°N. .Wind is 
the dominant mixing factor. The temperature is warm with a small season- 
al variation. The flow is generally zonal with a net westward transport. 
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Figure 19. Regions o£ the North Pacific Ocean [Tully, 1964]. 




Figure 20. Drift in the North Pacific Ocean [Tully. 1964]. 
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The Subtropic is the largest of the regions. The flow is zonal 
and relatively slow, as is sho'Aoa in Figure 20. The evaporation exceeds 
precipitation, causing an evaporation- driven convection mixing. The density 
structure is primarily a function of the temperature structure. From 
spring to autumn more heat is added during the day than is lost during 
cooling at night, so that in September the warmest surface water occurs. 
The region extends north to between Latitudes 40°N and 45°N. 

The Subarctic region is characterized by excess precipitation 
and as a result there is a distinctive salinity structure. Fresh water 
accumulates in the upper 100 meters with a halocline in the next 100 
meters where salinity increases about 1 percent. The density in this 
region is primarily a fxmction of salinity. Mixing is due to the v.’ind alone. 

The Pacific Arctic region also has an excess of precipitation 
and a salinity halocline. There is an ice cover in the winter and until mid- 
summer the heat added is used to melt the ice , followed by a rapid and 
shallow warming due to the long northern days. Convective mixing is 
dominant in this region. 

The Transport regions occur at the western boundaries of the 
ocean. An intensified meridianal flow transports water between the 
regions of zonal flow (Figure 20). The water does not remain in these 
transport regions long enough to fully acquire the regional characteristics. 

The Kuroshio Current is the western boundary for the Sub- 
tropic region. It transports a thick layer of warm saline water formed 
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in the Tropics near the Philippine Islands past the coast of Japan losing 
some of its heat along the way , but never becoming fully adjusted to the 
climate of the region. Its boundaries are easily identified due to large 
horizontal temperature gradients. 

The Oyashio Current is the western boundary for the Sub- 
arctic region. It has a smaller rate and volume transport than the 
Kuroshio Current. The flow is coast hugging and becomes strongest 
during the summer when it is swelled because of the ice melting in the 
Arctic. 

The Confluence Region is the area where the Kuroshio and 
Oyashio Currents meet. The resultant flow is eastward. The boundaries 
between these two flows can have temperature gradients as high as 5°C 
in a mile. Eddies from one stream can be dissolved in the other stream. 

As the Confluence moves east it looses its identity and is dissipated in 
the slow eastward trans-ocean flow (Figure 20). More than two years 
are required for the water to cross from Asia to America [Uda, 1963] 
so the surface processes can fully re-establish the regional characteristics. 

2. Microstructure Distribution 

The data available to this author necessitated that the regions 
investigated be restricted to the Subtropic , Subarctic , Kuroshio Current 
and Confluence regions. The first region and the one most thoroughly 
examined was the Subtropic for the purpose of comparing the results with 
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those o£ Roden [1971] who had conducted a detailed study o£ the vertical 
temperature and salinity structure in the central North Pacific in April 
1968 using STD data. 

Six locations were selected in the Subtropic region to see i£ 
the mean temperature perturbation vertical thickness scale was related 
within the region. The mean temperature perturbation vertical thickness 
scale was £ound £or each o£ the 68 available traces using a 1 meter Alter- 
ing interval. The average was then found for the traces from each of the 
six locations. This average mean temperature perturbation vertical 
thickness scale will be referred to as the average thickness scale. Also 
for each trace , the maximum amplitude of the temperature perturbation 
was found and then averaged for each location. Table II presents the 
results for the Subtropic region. The minimum, maximum, and average 
values for each location are listed along with the number of traces analyzed 
for each location. Roden's values are given for similar locations for 
comparison. The average thickness scale varies between 8.7 and 15.0 
meters which compares favorably with Roden's results which show a varia- 
tion from 13 to 16 meters for this region. One reason that a wide range 
of variations was found is that in analyzing the data , Roden used zero 
crossings to compute the thickness scale while in this analysis , inflection 
points from the digitized record were used because , as was pointed out in 
Section IV F, they represent layers resolved by the digitizing process. In 
addition , the 68 traces that were averaged represent all of the months of 
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the year and not just April. Another factor is that the maximum depth 
of the XBT traces was 750 meters while the STD record extended to 
1500 meters. 

The average maximum amplitude of the temperature perturba- 
tions varies between 0. 30 to 0. 49°C which is from the same minimum to 
three times the maximum value Roden found. This difference is probably 
largely due to the STD not recording the total magnitude of a temperature 
change since its response time is about one- third as fast as the XBT. 

The STD is also temperature compensated in the salinity sensor and large 
temperature changes in the water column cause a salinity spike that will 
show by comparison with the temperature record the whole change was 
not recorded. Another important factor ?ffecting Lhe STB is that it is 
not fi'ee from the vertical motion of the platform. This vertical motion 
as well as the physical size of the instrument could distort the temper- 
ature structure before it is sampled. 

Table III gives the results for the Subarctic region. Only 12 
traces were analyzed. The presentation is the same as for the Subtropic 
region with the exception of the omission of the Roden values since he did 
not have any results from this region. The average thickness scale varied 
between 17. 4 and 82 meters. The location for the maximum thickness 
scales were in the Gulf of Alaska in the Alaskian gyre. The three traces 
analyzed were for March, April, and July with a mean thickness scale of 
84m , 101m, and 51m respectively. The general behavior (more 
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microstructure activity present during the long summer days and less 
during March , the coldest month) o£ the conditions is as would be expected 
since the area is unstable and with the cold winter, convection mixing would 
tend to homogenize the water column and with the summer heating the 
microstructure would reform. The average maximum amplitude of the 
temperature perturbations varied between 0. 16 and 0. 42°C. 

The IQaroshio Current and Region of Confluence results are 
shown in Table IV. The only comment on the data is that at 40 °N 149°F 
the average thickness scale is smaller than for the rest of the region. 

This small value indicates a large concentration of microstructure activity 
and should be expected because this is the area where the Kuroshio and 
Oyashio currents meet , creating a region of intense vertical shear that is 
effective in producing microstructure. 

Table V gives results from two other areas. One location was 
in the Luzon Strait between Luzon and Taiwan. This also demonstrates 
that there is a large amount of microstructure present and again is in a 
location that current shear could be expected due to the restriction of 
flow. 

Also in Table V are the results of a 28 trace time series taken 
on two days in the California Current 200 nautical miles west of Monterey , 
California. Even though the extremes of the mean thickness scale are 
widely separated, 19 of the traces or 68 percent were within + 2 meters 
of the average value (13.0 meters). Table VI gives the local time, mean 
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thickness scale, and maximum amplitude o£ the temperature perturbations 
for each of the 28 traces. The hours are not consecutive because some of 
the traces were rejected by FNWC , Monterey \^hen received for process- 
ing. Without weather information, it is not possible to correlate the 
microstructure mean thickness scales and the local time. 

The average thickness scale data from Tables II, III, IV, and V 
has been plotted in Figure 21 with the latitude of the location, and shows 
quite clearly that in the north Pacific Ocean for an increase in latitude 
the average mean temperature perturbation vertical thickness scale in- 
creases. A similar plot. Figure 22, shows that the average maximum 
amplitude of the temperature perturbations decreases with an increase in 
latitude. Due to the scatter of data the relationship is not as clear as 
for the thickness scale. Figure 23 is a plot of the maximum amplitude 
of the temperature perturbations versus the average mean temperature 
perturbation vertical thickness scale. Although the scatter of data is 
wide, the relationship of the plot can be determined. For an increase in 
the thickness scale of the microstructui'e the maximum amplitude of the 
temperature perturbation will decrease. Sverdrup, Johnson, and Fleming 
[1942] have shown an increase in water instability \^^en preceding north 
from Hawaii to Alaska. With an increase in instability the thick layers 
of turbulent convective mixing and of small temperature gradient would 
become thicker, eliminating the thin layers of high gradient. This is a 
possible explanation for the results shown in Figure 21. 
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Climatological data was not obtained so that no conclusions 
could be drawn concerning the local time and the microstructure activity. 
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Subtropic Region 
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Subarctic Region 
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Transport Regions 




63 



TABLE V 



Other Areas 



California Current 


36-37N, 125-126W 


28 traces 


Mean Temperature 


MIN 


7.9 


Perturbation Vertical 


MAX 


19.3 


Thickness Scale (m) 


AVERAGE 


13.0 


Maximum Amplitude 


MIN 


0.16 


Of The Temperature 


MAX 


0.33 


Perturbations (*^0) 


AVERAGE 


0.24 


Luzon Strait 


20-24N, 120-124E 


2 traces 


Mean Temperature 


MIN 


6.5 


Perturbation Vertical 


MAX 


7.8 


Thickness Scale (m) 


AVERAGE 


7.15 


Maximum Amplitude 


MIN 


0.43 


Of The Temperature 


MAX 


1.00 


Perturbations (°C) 


AVERAGE 


0.71 



64 



TABLE VI 



California Current Time Series 



Local 

Time 


Mean Temperature 
Perturbation Vertical 
Thickness Scale (m) 


Maximum Amplitude 
Of The Temperature 
Perturbations (°C) 


11 Nov 68 0200 


15.0 


0.24 


0400 


15.4 


0.22 


0500 


12.8 


0.32 


0600 


13.3 


0.16 


0700 


16.6 


0.20 


0900 


13.1 


0.20 


1000 


15.9 


0.18 


1100 


15.0 


0.33 


1200 


10.7 


0.29 


1300 


11.0 


0.25 


13 Nov 68 2300 


7.9 


0.20 


14 Nov 68 0100 


12. 5 


0.22 


0200 


8.7 


0.16 


0400 


19.3 


0. 10 


0500 


13.5 


0.16 


0600 


15.7 


0.18 


0800 


10.7 


0.22 


0900 


13.1 


0.20 


1000 


9.9 


0. 24 


1100 


13. 5 


0.16 


1200 


12.9 


0.14 


1300 


11.0 


0.20 


1400 


10. 5 


0.27 


1500 


13. 3 


0.20 


1700 


13.1 


0.22 


1800 


14. 3 


0.22 


1900 


12.8 


0. 26 


2100 


12.4 


0.24 
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Figure 21. Plot o£ the Average Mean Temperature Perturbation 
Vertical Thickness Scale versus Latitude showing an increase in 
the vertical thickness scale for an increase in latitude between 
20°N and 55°N. ( ^ , points along Longitude 156-157^W; • , all other 

points). 
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Figure 22. Plot of the Average Maximum Amplitude of the 
Temperature Perturbations versus Latitude showing the Amplitude 
decreases with an increase in latitude between 20°N and 55°N. 

( A , points along Longitude 156-157°W; • , all other points). 
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ature Perturbation Vertical Thickness Scale increases between 
20^N and 55^N. ( A , points along Longitude 156-157^W; • , all other 

points). 



68 



V. OPERATIONAL ANALYSIS 



A. OPERATIONAL REQUIREMENT 

Underwater acoustic energy transmissions are extremely vital to the 
Navy , especially in the anti-submarine warfare (ASW) field. Some of the 
military uses and some of the non-military uses of underwater sound as 
reported by the National Academy of Sciences [1970] are listed below by 
frequency bands. 

1. Less than 100 Hz 

a. Passive and active long range sonar 

b. SOFAR navigation and positioniiig 

c. Deep exploration and mapping 

2. 100 to 5000 Hz 

a. Shallow seismic exploration and mapping 

b. Warning bells 

c. Scuba diver communications 

d. Passive and active ship sonars 

e. Mines 

3. 5 to 15 kHz 

a. Underwater telephone transmissions 

b. Depth sounding equipment 

c. Fish finding 

d. High frequency ship sonar 
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4. 15 to 100 kHz 



a. Short range underwater telephone transmissions 

b. Depth soundings 

c. Fish finding 

d. Ship positioning 

e. Object location 

£. Low resolution scanning sonars 

g. High resolution sound velocity profiHng 

h. Homing torpedoes 

5. Greater than 100 kHz 

a. Doppler sonar for navigation or docking 

b. High resolution scanning sonar 

For low frequency acoustic transmissions thermal microstructure 
should not have a significant effect on the propagation because the acoustic 
wave length is large relative to the effective radius of a patch of thermal 
inhomogeneities. However, losses at low frequencies have been related to 
scattering by the index of refraction inhomogeneities [Urick, 1967]. Above 
a few hundred cycles the microstructure becomes increasingly important 
and affects sonar transmissions adversely. 

Sagar [1960] , by taking temperature profiles while collecting acoustic 
(14. 5 kHz) transmission data, was able to associate acoustic intensity 
fluctuations with the presence of micro structure while in the absence of 
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microstructure the variation was small and independent o£ the Mintzer 
dependence on range. This illustrates the necessity for the Navy to be 
concerned with thermal microstructure. 

B. OPERATIONAL METHODS. 

XBT data is being used by the Navy to obtain the gross temperature 
profile for use in the solution to ray path equations and for diagrams to 
predict sonar performance. This can be accomplished on board a ship 
or can be accomplished by Fleet Numerical Weather Central , Monterey 
in a more sophisticated manner, and the products supplied as a service to 
fleet units. 

On board ships the parameters most used are the sea surface temper- 
ature and layer depth. Tables can then be utilized to find the predicted 
range for active and/or passive sonar contacts. Little information con- 
cerning the temperature structure is presently being used in weapon 
selection or in determining the expected performance of the weapon 
once in the environment. 

FNWC , Monterey supplies information on request to support ASW 
operations. The products include sonar range predictions and propagation 
loss predictions. The major factors vdiich are used for the prediction 
computations are water mass type , bottom type , and depth of the water. 
Temperature information is furnished from the historical data bank with 
no real time input from the area in question. Historical data is used to 
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prevent the computer solution from being biased in case the ship XBT was 
from a small patch of water rather than being a representative sample of 
the total area. Again, no microstructure information is being used. 

C. POSSIBLE UTILIZATION OF MICROSTRUCTURE INFORMATION 

ON AN XBT TRACE. 

A ship has the capability of sampling the temperature in its area wlien 
desired. However , the XBT trace is not being fully utilized. By closer 
examination of the trace , the sonar operator could derive an indication 
of the quantity of microstructure present in his operating area. With this 
additional knowledge the operator would know whether to expect a variation 
in sonar reception or to expect steady reception. He would also know of 
the possibility of picking up contacts beyond the predicted range that 
might be held only intermittently if microstructure is present. He could 
request reduced range to maintain contact if the contact is gained at the 
predicted range and microstructure is present. Use could also be made 
of microstructure information in selecting the range for launching a homing 
torpedo. The variation in acoustic intensity will be larger for the frequen- 
cy range of the homing to3rpedo than it will be for sonar frequencies , so 
it would enhance the opportunity for success of the torpedo if the presence 
or absence of thermal microstructure is considered. 

FNWC , Monterey could aid the ASW problem by predicting the relative 
magnitude of microstructure to be expected for an area by examining the 
microstructure on digitized traces now in the data bank for the area. This 
could be made a part of the sonar range prediction now ava ilable. 
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VI. CONCLUSION 



A. SUMMARY 

The motive behind the investigation o£ the thermal microstructure 
distribution in the north Pacific Ocean pursued in this thesis has been to 
determine -what information can be obtained from XBT records and then 
used to give some prediction of the influence on sonar transmissions and 
receptions. Various techniques were explored to determine a measure of 
the thermal microstructure on an XBT trace ending with the selection of 
a modified binomial filter, suggested by Roden [1971], being preferred. 
Data was extracted from. Fleet Numerical Weather Central, Monterey's 
XBT data bank, for four of Tully's [1964] climatic regions of the north 
Pacific Ocean. 

B. CONCLUSIONS 

Based on the examination of the data from the Subtropic, Subarctic, 
Kuroshio Current , and Confluence regions , the following conclusions are 
offered: 

1. The modified binomial filter is a technique that will produce a 
mean temperature perturbation vertical thickness scale for the thermal 
microstructure found on an XBT trace. 

2. The digitizing process used by Fleet Numerical Weather Central, 
Monterey, divides the trace into layers of linear tempera ''rure gradient 
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so that the use of the digitized inflection points will give the same result 
as that obtained using inflection points from the temperature perturbation 
curve. 

3. The binomial filter will also produce a measure of the amplitude 
of the temperature perturbations. 

4. The average mean temperature perturbation vertical thickness 
scale increases as latitude increases, vlaile the average maximum ampli- 
tude of the temperature perturbations decreases as latitude increases 
in the north Pacific Ocean between 20N and 55N. 

5. The average maximum amplitude of the temperature perturbations 
decreases if the average mean thickness scale is increased in the north 
Pacific Ocean. 

6. The average mean temperature perturbation vertical thickness 
scale was found to be more dependent on the latitude than TuUy's Climatic 
Regions. 

7. No association could be determined between the time of day 
or season and the mean thickness scale. 

C. RECOMMENDATIONS 

The XBT digitized data could be used for a variety of studies. Recom- 
mendations for further study enumerated below are specifically related to 
the development of the data processed in this thesis. 

1. One location should be selected for an extensive study to deter- 
mine if the XBT data will reveal a correlation with time of day and season. 
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Two likely areas in the Pacific Ocean for such a study are the San Diego 
and Hawaiian operating areas used by the Navy. XBT data should be very 
plentiful from these areas. 

2. A study should be conducted to determine the sonar transmission 
and reception variation that can be expected from the influence of thermal 
microstructure that appears on XBT traces taken during the study. 

3. Fleet Numerical Weather Central , Monterey , California , should 
investigate the feasibility of including a microstructure quantitative pre- 
diction as a part of the Sonar Performance forecast for the active sonars. 

4. Fleet Numerical Weather Central , Monterey , California , should 
add to the identification portion of each digitized XBT trace the climatol- 
ogical data ^ such as wind speed, sea surface temperature, and cloud cover, 
for ease in evaluation of the trace. 
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